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ABSTRACT 


Physical and thermodynamic concepts are used to develop a potential function 
for application to high temperature polycrystalline material response. Inherent 
in the formulation is a differential relationship between the potential function 
and constitutive equation in terms of the state variables. Integration of the 
differential relationship produces a state variable evolution equation that 
requires specification of the initial value of the state variable and its time 
derivative. Analysis of these conditions showed that the initial loading rate, 
which is directly related to the initial hardening rate, can significantly 
influence subsequent material response. This effect is consistent with observed 
material behavior on the macroscopic and microscopic levels, and may explain the 
wide scatter in response often found in creep testing. The material used for 
the study, cast and wrought IN100 at 732°C ; was tested in tension at different 
strain rates, creep, stress relaxation, and reversed inelastic flow. 

INTRODUCTION 


The research presented in this paper is directed toward developing a state 
variable constitutive model for metals in high temperature environment where 
rate effects are important. The underlying concept used in the model is to 
develop a consistent system of equations to predict the inelastic strain rate 
and evolution of the state variables that are derivable from a potential 
function. The essential structure of the theory is based on the maximum 
plastic work inequality, the rate of work hardening inequality, and dislocation 
dynamics . 

The maximum plastic work inequality, expressed using the concept of a yield 


surface, can be written as (refs. 1 and 2) 
(o ij ' £ °- 


d.i) 
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The quantity de^ is the inelastic strain increment due to the stress step 

a , - o° . The stress a., is on the current yield surface and a° . is any other 
ij ij ij ij • 

stress either inside or on the yield surface. The inequality requires that the 
inelastic strain increment is normal to the yield surface and that all yield 
surfaces are convex. Eq. (1.1) is an axiom resulting from observed plastic 
behavior of metals. It is not derivable from thermodynamics and is not 
necessarily valid for all materials. It does, however, capture the essence of 
how many metals behave. 

The property of work hardening, as identified with a large class of metals, 
can be further characterized. In these materials it has been consistently 
observed that the stored energy of cold work, W^, increases with deformation at 
a decreasing rate during a continuous deformation history. This can be expressed 
by 

dw!) - dwj <_ 0 , (1.2) 

li h 

where dW and dW_ are the increments of stored energy associated with identical 

^ Yi h 

strain (or stress) increments in two different stored worked states W^ and 
such that (ref. 3). The exact structure of a cold worked state is not 

fully understood, however, it is generally«accepted that the energy is stored 
through the development of a system of dislocations (ref. 4 and 5). 

In the initial state of deformation dislocations multiply and tend to 
arrange themselves into groups or clusters occupying only a small fraction of 
the material volume. Most of the deformed material is dislocation free. As 
deformation increases the dislocation clusters form continuous walls separating 
relatively perfect cells. With further strain, secondary slip systems are 
activated and the density of dislocations increases with other point defects 
appearing. The cell size decreases at a decreasing rate as strain increases 
and tends to stabilize. The cell size and dislocation density at any time are 
influenced by the initial microstructure, temperature and loading history up 
to that time. In addition, it is important to recognize that a substantial 
portion of the character of the microstructure is established early in the load 
history when the observed macroscopic inelastic strain is very small. 

In some cases, the stored energy is partially recovered in time at elevated 
temperature or fully recovered through recrystalization. This effect could be 
important for metals at high, temperature that have been work hardened in the 
initial configuration. This situation is typical of cast and wrought super- 
alloys. In the case of recovery, W^ < W^, and the stored energy increments dW^ 
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and dW^ are both negative. In this case, eq. (.1.2) requires the rate of 
change of recovery to be positive. That is, softening occurs at an increasing 
rate; conversely, hardening occurs at a decreasing rate, and both hardening 
and softening tend to stabilize. 

The maximum plastic work inequality and the rate of work hardening 
inequality satisfy most observed response in metals at high temperature under 
continuous load histories. However, the essence of the theory can be easily 
extended by a broader class of materials and load histories by requiring the 
two inequalities hold together rather than individually; that is 

(a.. - cr?.)de*. - (dw!) - dwj) > 0 . (1.3) 

A similar relation was suggested by Ponter (refs. 5 and 6) and the consequences 
are extensive. First, since the stored energy is generally small compared to 
the total plastic work during a typical deformation history, the restriction 
on the hardening or softening rate (eq. (1.2)) can be softened in some cases. 
This allows, for example, a jump in the rate plastic working to produce a jump 
in the rate of work hardening that is not restricted by eq. (1.2), The effect 
appears possible in situations where there are jumps in strain rate early in the 
development of microstructure. 

Another essential feature of the constitutive model is to use the result 

of Rice (ref. 7) showing that the components of an inelastic strain rate tensor, 

• I 

e^, are derivable, at each instant during the deformation history, from a- 
potential function, of the stress, i.e. 

.1 = aflfg, historyl 

ij do.. ' (1,4) 

This result is based on the physical notion of conventional crystalline deform- 
ation: At a given slipped state, the rate of permanent shearing on a particular 

slip system depends on the stress at that point only through the shear stress 
acting on the slip system. Thus, for a given prior history, eq. (1.4) was 
shown to give time dependent stress strain equations consistent with both 
plausible macroscopic and microscopic idealizations. More recently, Ponter 
and Leckie (ref. 8) extended the formulation to polycrystalline metals at high 
temperature. Following the methods in ref. 6 they established a potential 
function for a constitutive equation of the Baily (ref. 9) and Orowan (ref. 10) 
type that contains one state variable, s, to describe the hardness or stored 
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energy of the material. In the Ponter-Leckie development the assumption of a 
local (microscopic) potential was essential to establish a macroscopic potential 
that can be used to derive both the inelastic flow equation and state variable 
evolution equation. The Baily— Orowan the macroscopic flow equation contains 
an over stress argument of the type - s, where $ is a scalar function of the 
stress tensor. 

Another state variable constitutive equation that uses a single state 
variable to describe the hardness or resistance to inelastic flow was developed 
by Bodner and Partom (refs. 11 and 12). Inherent in the representation is the 
absence of a yield surface. That is, the inelastic strain rate is non-zero 
for all non-zero values of stress. This is a continuous flow equation without 
separate loading and unloading conditions. It has been used to successfully 
predict the response of two superalloys at high temperature (refs. 13 and 14). 

In part, the success may result from using the entire load history rather than 
excluding the load history inside the yield surface. For example, recall that 
the character of the dislocation substructure is established very early in the 
deformation history prior to classical yield. 

Another key feature of the Bodner-Partom model is that the specific form 
of the flow equation was motivated by dislocation dynamics. Both of the well 
accepted representations for dislocation velocity, developed by Gilman (ref. 15) 
and Vreeland (ref. 16) are embodied in the Bodner-Partom formulation. One 
shortcoming of the Bodner-Partom representation is that the basic structure 
of the evolution equation for the hardness state variable is developed by 
phenomenological methods. The evolution equation is consistent with the observed 
properties of stored energy (eq. (1.2)) but does not possess a formal 
mathematical or physical derivation. 

Specifically, in this paper it is shown that the concept of a macroscopic 
inelastic potential function is compatible with the essential features of high 
temperature material response expressed in eq. (1.3) . A potential function 
concept is then used to derive a state variable evolution equation directly from 
the inelastic flow equation. A specific example is developed using the Bodner- 
Partom representation. It is shown that the Bodner-Partom evolution equation 
has a mathematical structure very similar to the derived evolution equation. 
However, the derived evolution equation has a new important property. The 
derived representation depends on the initial hardening rate that can vary from 
test to test for the same initial material microstructure. The representation 
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is used to predict the response of IN100 at 732°C. The material constants are 
determined from tensile response data and the experimental response in, creep, 
stress relaxation and reversed inelastic flow is analyzed. 

II DEVELOPMENT OF A POTENTIAL FUNCTION FORMULATION 

Let us introduce a state variahle Z as a macroscopic measure of the effect 
of the dislocation microstructure on deformation. It is designated as the 
hardness or resistance to inelastic flow and the units are that of stress. 
Similarly a macroscopic strain state variable, e* 1 , is defined as a measure change 
in geometry associated with, the development of dislocation microstructure. 

The variables Z and e* 1 are scalar functions of tensor valued arguments and 
defined so that an increment of stored energy of cold work can be calculated as 

dW* 1 = Zde h = de h dt (2.1) 

It is required that the stored energy increment, dW* 1 , and hence the valuables 
Z and e , satisfy the conditions outlined in the Introduction. The increment of 
inelastic work associated with an increment of inelastic deformation is 

dwI - V e u " Vu dt (2 - 2) 

In a constant hardness state the maximum plastic work inequality (eq. (1-1)) 
can be written as 

• I 

E ijd a ij — Z = const. (2.3) 

Similarly, the rate of work hardening inequality, for constant stress eq. (1.2.) 
can be rewritten as 

• h 

e dZ 0, = const. (2.4) 

and the combined inequality (eq. (1.3)) is 

e. . da. . - e h dZ > 0. (2.5) 

ij ij — 

Z = const. a.. = const. 

Using the result of Rice (eq. (1.4)) there exists a potential function 4>(o,Z) 
such that 

•I 9$ (a, Z) 

G ij 3o“ * (2*6) 
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Further, using the eqs. (2.5) and (2.6) , and the result of Ponter (ref. 17), 
let us assume that for <j> to exist and be physically acceptable, that 


•h 3$ Co, ZD 
3Z 


(2.7) 


and 


3ifr 

3a 


ij 


da ij + If dz - 


> 0. 


Therefore, <p exists and it is convex in the space (a, 2} (ref. 16). The existence 
of the potential function <f>(a, Z) is consistent with the essential structure of 
work hardening or softening metal plasticity. As a result of the existance of 
<f> and from eqs. (2.6) and (2.7) it follows that 



3a. . 
ij 


3et. 

3Z 


( 2 . 8 ) 


Using the idea of an instantaneous tangent modulus in one dimensional 
plasticity, the stress rate state variable Z can be related to the strain 
rate variable e by the function h; i.e. 


he 


" h 3Z * 


(2.9) 


That is, h is the instantaneous slope of the Z ~ e response curve. For the 

following analysis let h = h(Z,...) be independent of the current stress. 

Finally, eq. (2.8) can be used to derive the state variable evolution 

• I • I 

equation directly from the flow equation. Let us assume e = e „ (a, Z) is 
specified, then from eq. (2.8) and (2.9) 

dZ = - h — u - 


3Z 


da ij + df (Z,...) 


( 2 . 10 ) 


where f is independent of stress. Integrating on the time interval C0,t3, the 
evolution equation for the hardness can be written as 
. a (t) Se*. 

Z (t) - Z Q =-h/ a 1J (0) da + Cf(Z,...) - f'(Z 0 ,...)3 . (2.11) 

The integral can be evaluated for any specific choice flow law involving only 
one state variable. 
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Embodied in eq. (2.11) is an important property that is not widely 

• • 

recognized. The parameters Z^ = Z(0) and = Z(0) are initial conditions. The 
parameter Z^ describes the initial state of the microstructure that, for example, 
would be the same for all specimens from the same heat of material. However, Z^ 

is the initial rate of hardening of the microstrucutre and would be expected to 
vary from test to test, depending on the initial loading conditions. 


Ill Bodner-Partom Constitution Equations 


The Bodner-Partom constitutive equation (refs. 11, 12, 18) is a fully 
developed three dimensional theory that has even been extended to anisotropic 
materials (refs. 19 and 20). Since the main purpose of the paper is to evaluate 
the properties of the evolution equation, it is convenient to restrict the 
analysis to a one dimensional fora of the constitutive equation. This is 
consistent with the accompanying experimental program. 

The isotropic constitutive theory of Bodner and Partom is based on the 
assumption that the total strain rate, e(t), can be separated into elastic, 

• 0 t J 

e (t), and inelastic, e (t), components. Let E represent the elastic modulus, 
then the Bodner-Partom equation can be written in a one dimensional fora as 


eCt) 


g(t) 

E 


+ e X (t) 


(3.1) 


where o(t) is the current value of the stress. Inherent to the theory is that 
the inelastic strain rate is non-zero for all non-zero values of stress. The 
specific representation used by Bodner and his co-workers for the inelastic 
strain rate is given by 


•I 

£ 


(t) 


= D 


o(t) 

0|cr(t) 


r (n+1) /Z^n-, 

esp c- — <— ) : 

c 


(3.2) 


The constant D~ represents a limiting value of the inelastic strain rate and is 
generally taken at 2x10 / S3 sec unless the strain rates are very high. The 
constant n controls the strain rate sensitivity and Z(t) is the hardness state 
variable. The general mathematical structure of eq. (3.2) is based on dis- 
location dynamics expressed in the context of continuum mechanics and has proven 
consistent with the observed response for many metals. The formulation is 
similar to the classical yield surface theory. The structure of the Prandtl-Reuss 
formulation is preserved, but a yield surface itself is not part of eqs. (3.1) 
and (3.2). 
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Necessary for the use of the above equations is a representation for the 
state variable Z. The specific representation proposed by Bodner et. al. 

(ref. 18) and used for superalloys (refs. 13, 14 and 21) is based on the concept 

•I 

that only the inelastic rate of working W and the current hardness control the 


rate of hardening. The representation is written as 

•P Z_Z 2 r 

Z = m(Z - Z)\t - AZ. (.-—■) 

1 1 \ 


(3.3) 


with Zq designated as the inital value of Z. The two terms in eq. (3.3) are 
defined so that AZ^C (Z-Z^J/Z ^3 is negligible during rapid load histories. 

Thus, during a tensile test that is fast compared to creep test, eq. (3.3) 
reduces to the first term alone. The constant Z^ corresponds to the maximum 
value for Z and m is an exponential coefficient controlling the rate of harden- 
ing. For a long time response, such as creep, a second term corresponding to 

•P 

hardening recovery is necessary. During the minimum creep rate both e and a 

are constant, thus Z is constant (Z = 0) and the rate of hardening must equal 

the rate of recovery. The coefficient Z^ corresponds to the minimum recoverable 
value of hardness, and A and r are the coefficient and exponent, respectively, 
controlling the rate of hardening recovery. 

At this point in the paper it is now possible to carry out one specific 
objective: Derive an evolution equation for Z by using eq. (3.2) in 

eq. (2.11) and to compare the result to eq. (3.3). To begin, let us define 

l/2n 

(3.4) 


B “ D 0 <£r> 


, D n+1 /Z.2n 
and R = -r — (-) 
zn a 


so that the integral in 
I = B / R- 1/2n e~ R 


eq. (2.11) 
dR . 


can be written as 


(3.5) 


For continuous histories eq. (3.5) can be integrated by parts N times to 
obtain a series representation. Letting p = l/2n, the series can be written as 


I(R) = B e _R R P C- 


R~ 


(1-p) (1-p) (2-p) T (1-p) (2-p) (3-p) 


+ ... 3 


(3.6) 


which converges for all R on (0,«). The function I(R) is almost constant for 
large values of R (small values of stress) and decays with decreasing R (increas- 
ing stress). In this application of eq. (3.6) the constants D Q and n are 
generally known so eq. (3.6) is totally defined. 

Finally, since eq. (3.5) could be integrated to obtain eq. (3.6), then 
eq. (3.6) is an exact differential. This implies that the contribution of 
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I(R) to the current hardening rate depends only on the current value of the 
state, Z, and the current stress. Thus, for all histories producing the same 
state and stress the integral term will have the same contribution to the 
hardening rate. 


IV PROPERTIES OF THE TWO EVOLUTION EQUATIONS 

The derived evolution Equations, (eqs. (2.11) and (3.6)), can be put in a 
form similar to the Bodner— Partom eq. (3.3). Noting from eqs. (3.2) and (3.4) 
that I(R) can be written as 

(4.1) 


I(R) = g(f)e Z 


where 


✓£>. _ n /n+lxl/2n -p r R R 2 , 

S V “ D 0 ( 2n J R L (l-p) + (1-p) (2-p) + (1-p) (2-p) (3-p) 


+ ...3 (4.2) 


and R(^0 is given by eq. (3.4). Thus, the derived evolution Equation, 
(eq. (2.11)), becomes 

z = Z 0 - hCg(f) - ggCf)^ 1 + Cf(z,...) - f(z 0 , ...):. 


(4.3) 


The Bodner-Partom Evolution Equation (eq. (3.3)) can be rewritten for convenience 
as 

Z-Z 

Z = m(Z 1 - Z)ae I - AZ^-^)* . (4.4) 


Direct comparison of eqs. (4.3) and (4.4) shows a very strong mathematic 

similarity. Both equations have two terms of similar type. One term is linear 

in the inelastic strain rate, e Z , and the other term is independent of the stress 

a. The derived evolution equation has the initial hardening rate term Z^. The 

coefficient of the strain rate term in the derived equation is a product of a 

2 

modulus h and a non-linear function of (— ) ; whereas, the corresponding term in 
the Bodner-Partom is linear in both a and Z. The function f(Z,...) in the 
derived equation is consistent with logarithmic function of Z in the Bodner- 
Partom equation. In general, eq. (4.4) can be considered almost a special case 
of eq. (4.3). 

To evaluate the properties of the derived evolution equation further it is 
convenient to develop a representation for the function f(Z,...). A power series 
expansion of f in Z would produce a series of exponential time terms when 
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eq. (4.3) is solved. Since the general solution is exponential in nature, a 
reasonable approximation is to use a single exponential term. Thus, let us 
assume a linear representation, namely 

f(Z) - f(Z Q ) = -o(Z - Z Q ) (4.5) 

where a is a parameter that is independent of Z and a. Further, in order for 
the representation for Z to approach a stable asymptotic value it is necessary 
for a > 0. It is expected that this representation should contain many of the 
essential features. 

Even though the mathematical structure of eqs. (4.3) and (4.4) is similar, 

the physical interpretation of the terms is different. First notice that t = 0, 

• • • 9 

Z = 0 in eq. (4.4) whereas Z = Zq in eq. (4.3). It is expected that the initial 

hardening rate should depend on the initial loading conditions, that is 

• • • • • « • 

Z 0 = Zq(£q) = Zq(Oq/E) where and 0 q are the initial strain and stress rates. 
This adds considerable flexability to the model and is consistent with the 
initial formation of dislocation microstructure as described earlier. Next, the 
strain rate term in the third term of eq. (4.3) produce basic tensile response 
properties and the stable value for Z. The strain rate term of eq. (4.3) and 
the second term of eq. (4.4) characterize the long time recovery response 
properties of the model. 

V EXPERIMENTAL PROGRAM 

Sixteen mechanical tests have been conducted on IN100 at 732°C (1350°F) at 
the Air Force Wright Aeronautical Laboratory, Ohio and Mar-Test Inc., Cincinnati, 
Ohio. The material was obtained at different times from different heats resulting 
in five groups of specimens designated as series C, G, T, GT and ENTEN. The 
experimental program, summarized in Table 1, includes eight tensile tests, seven 
creep tests and one combined test. The controlled experimental variable is 
shown in Table 1 and the observed stable values for stress or secondary creep 
rate is also given for the tensile and creep tests, respectively. 

The results of seven tests conducted under constant strain rate control 

_3 

ranging from 1.4x10 sec to 1.6x10 sec and one test under constant head 

-4 -1 

rate control at 8.3x10 sec are shown in fig. 2. There is significant 
variation in the level of the stress response due to the imposed variation in 
strain and head rate. Note, however, the total accumulated strain in these tests 
is not important since several of the specimens were not failed. For four 
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different values of strain rate (Tests 2, 4, 6 and 8) the response obtained a 
maximum stress and maintained that value of stress for all subsequent values of 
strain. However for Tests 5, 7 and 9 a different response was obtained. In 
these experiments the value of the stress decayed from the maximum value obtained 
at about one percent strain. In test 5 the amount of reduction in stress to a 
lower stable value was small. In Test 7, at a lower strain rate, the reduction 
in stress to a stable value was larger; and in Test 9, the stress did not 
stabilize at a lower level. This wide variation in response might arise since 
the eight specimens are from four different heats. Since both types of 
response were observed at both AFWAL and Mar-Test Inc., it cannot be accepted 
as an experimental problem. 

The results of six creep tests are shown collectively in fig. 6b for 
up to 100 minutes. The variation in the creep stress was almost twofold 496 
to 896 MPa (72 td 130 KSI), and the corresponding minimum creep rates are given 
in Table 1. Although not shown. Test 10, 14 and 15 obtained tertiary creep at 
times ranging from 100 minutes to about 1200 minutes. The response curves do 
not exhibit a significant primary creep phase and most of the response is 
dominated by tertiary creep. This is typical of other superalloys. Test 20, 
creep at 896 MPa (130 KSI) shown in fig. 6b, is not ordered with respect to 
the other tests; however, this could result from speciments variation or loading 
conditions as discussed later. 

If the deformation mechanism controlling the tensile tests and creep tests 
are the same, then the controlled and observed variables from the two types of 
tests shown in Table 1 correspond to the same deformation process. That is, the 
stable value of stress obtained in a strain rate controlled tensile tests should 
correspond to the creep stress with the same constant (secondary) creep rate as 
the tensile test. A plot of the observed and controlled variables for both the 
creep at tensile tests is given in fig. 1. Considering the data is over five 
decades of strain rate, there appears to be reasonable consistency between both 
types of tests. Tensile Tests 7 and 8 correspond very closely to creep Tests 19 
and 18, respectively, as shown in fig. 1 and table 1. Thus, it does appear that 
the same basic deformation mechanisms control both creep and tensile behavior 
between 482 and 1100 MPa (70 and 164 KSI) at 732°C (1350°F). 
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VI APPLICATION TO IN100 

The next step in predicting the response of IN100 at 732°C is to determine 
the remaining material parameters in the evolution equation. This exercise is 
divided into two steps. First, constant strain rate histories are analyzed and 
then the results are extended to piecewise constant strain rate histories. 

Constant strain rate histories. The evolution equation that is derived 
from the potential function can be written as 


z = z Q - ct(z-z Q ) + hCi(o)-Kt): 


( 6 . 1 ) 


for continous deformation histories using eqs. (4.5) and (4.1) in eq. (4.2). 

During each constant strain rate tensile test the last term approaches a step 

to the steady value, A . An approximation of the hehavior can be made by 

s s 

replacing the integral term in eq. (.6.1) by a step function. Thus, assume 
eq. (6.1) can be approximated by 


Z = Z n - a (Z— Z ) + A pCt-t 3 
0 0 ss n 


( 6 . 2 ) 


where t is the time of the unit step pCt-t 3. The solution of eq. (6.2) for 
n n 

the hardness variable yields 


Z„+a Z n 

z - z o e- at + C-i—O) a-e-“ t > + 


ss r- ^ -t f, -a ( t- 1 ) v 

yLt-t 3 (1-e n ) 

a n 


(6.3) 


which is expected to be a reasonable approximation for times up to t R and at 
times much larger than t . 

Since a significant change in response occurs early in the history, let us 

define t as the time to reach the ultimate stress under a hypothetical elastic 
n 

loading conditions using the initial strain rate; ie. 


"n Ee, 


(6.4) 


where and are the ultimate stress and constant strain rate during each 

tensile test. The elastic modulus is denoted by E. Thus, the response up to 

time t characterizes the development of the microstructure before the onset of 

significant macroscopic inelastic deformation. The value of hardness, Z , 

n 

expressed from eq. (6.3) as 


z o + IT (1 - e " atn) 


(6.5) 
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( 6 . 6 ) 


and the steady state value (Long time) value is 
Z 


Z. + aZ. + A 
0 0 ss 


ss 


a 


Equations (6.5) and (6.6) can be used to evaluate the parameters in the evolution 

equation at three points: the initial condition Z^ and Z^; the onset of 

significant inelastic flow, Z ; and in the fully saturated condition, Z 

n ss 

The values of Z and Z can be calculated directly from the tensile data 
n ss 

for each test. Inverting eq. (3.2) and evaluating at time t R gives 


r -2n o / 0\-il/2n 

z - * "u'sh: ln( ^ ):i 


n 


(6.7) 


likewise at saturation 


Z = a £ n (_0):] 1/2n 

ss ss n+1 Dq 


( 6 . 8 ) 


Comparing eq. (6.1) to eq. (6.2) and using eq. (6.8) yields 

Z 


A = hCl(O) - I(— )1 
ss o 

ss 


which can be used with eq. (6.6) to define the parameter h as 


h = 


“< Z ss- Z 0>- Z 0 


:i(o)-u^): 

ss 


(6.9) 


( 6 . 10 ) 


The parameter h, as, expected, is a function of history noting from eq. (6.8) 

that (Z /c ) is dependent only on the imposed constant strain rate £». The 
s s s s u 

initial hardening rate can also be determined from eq. (6.5) 


z o- 


a < z „- z 0> 

1-e at n 


( 6 . 11 ) 


noting that Z can be determined from eq. (6.7). 


n 


The time parameter a can by estimated from the time required to obtain the 

— ct t 

steady state response. Assuming e vanishes at at = 5 and using saturation 
strain e^, then ae^/leQl =5. In this example 

a - C 0 |e 0 | , (6.12) 
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where = 10 sec 4 and the time response is also dependent upon the applied 
strain rate. The remaining parameter Z^ is choosen by interactive computation 
to obtain the best fit for the lowest strain rate response curve. The permits 
including the recovery effect as well as strain hardening. This, analysis 
produced = 3275MPa(475 KS1) . The representations derived directly from the 
tensile data for use with eq. (6.7) and (6.8) are 


o^ - 6.8948(218.79 x “ 0 - 77934 _ 58 .65 x -2 * 9144 ) 


_ _ r QQ/oricn n/ Q -0.31697 9.6381 N 

a = 6.8948(162.048 x - 0.047656 x ) 


(6.13) 


where x = log e/log and stress is in MPa(KSI). The representations were 
established to account for two different deformation mechanisms above and below 
approximately 900 MPa (see fig. 6). The other constants for IN100 used in the 
study were taken from ref. 14. They are 

E = 1.50x10"* MPa, n = 0.70 and = 2xl0 4 //3 sec 1 . 

Piecewise Constant Histories. The representation for constant strain rate 

histories can be summarized by substituting eq. 6.10 into eq. 6.2 and identifying 

Z /a from eq. (6.8) as a function of the applied strain rate, e_, to get 
s s s s u 

•7 

Z = Z - a(Z-Z ) + Ca (Z -ZJ-Z 3 . (6.14) 

o o ss o o cf(oT-uf^i 

• • 

Let us consider a step change in strain rate from to at time t^. 

From eq. (3.1) it is seen that this produces a jump in the stress rate, 

Aa 1 = EC(e 1 -e 0 )-(e^-e^)3 (6.15) 

but not a jump in the stress. Further it is expected that the inelastic strain 
rate would remain the same; therefore, from eq. (3.2), the value of the hardness 
Z is the same immediately before and after the jump in the strain rate. This 
implies that the microstructure does not change instantaneously. Note that 

eq. (6.15) also shows that separate unloading conditions are not required. 

• • • • 

If = ~ € q> then Aa^ = -2 EEq, which is very large compared to a stress rate 
near zero before the jump if the material is in the plastic range. Also, note 
that since the one test with reverse plastic flow does not indicate a significant 
amount of Bauschinger Effect, it is reasonable to assume an "isotropic hardening 
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rule" and use the same hardening rate equation for loading in both tension and 
compression. 

It is expected, however, that a jump in strain rate would produce a jump 

• . . 
in the hardening rate, Z. From eq. (6.14), a jump in the value of e^ to e^ 

• • 
should produce a jump in I(e Q ), z gs and “• It is necessary to update I(Eq) to 

1(e) to maintain the asympotic value for Z. In this case, as time becomes 

X * 

large, Z approaches zero for any choice of Z and a. Thus, Z and a define 

s s s s 

the new microstructure corresponding new strain rate and time required to 
arrive at the new structure. For IN100,it was found that updating Z gs , 
especially after several strain rate steps, produces erroneous results and 
estimating Z^ from the initial strain rate produced much better results. 

Thus, due to the lack of more specific information, Z gs was determined from 
the initial value of strain rate. This implies that the microstructure 
generated during the initial load secquence establishes some of the major 
features of the microstructure for all subsequent deformations. This appears 
to be reasonable for the simple deformation histories used in current study; 
but, may not be valid for all materials and all deformations. Finally, from 
the lack of better data during creep and stress relaxation testing, as 
fully explained later, it was decided to use the same value for a the entire 
history as a first approximation. Thus, for step changes in strain rate, the 
hardening rate described in eq. (6.14) is used with the values of a and Z 

ss 

determined by the initial strain rate. 

Extension to continuously varying strain rate histories could perhaps 
be modelled by using the current value for strain rate as described above. 

However, in the current paper only piecewise constant strain rate histories 
were used in the experimental program so it is difficult to test this 
hypothesis. Extension to other histories and materials requires additional 
development along with a better experimental description of micros true tural 
changes . 

VII CALCULATED AND PREDICTED RESULTS 

The observed and calculated tensile response curves are shown in fig. 2. 

The calculated curves match the observed response with remarkable detail. The 
constants were determined with only two parameters, Z^ and Cq, being determined by 
computational techniques. This was done using Curve 9 only. Further only 
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tensile data was used to determine the constants in eq. (6.13). All other 
predictions were made using only these constants. 

It is interesting to note that the shape of both the observed and calculated 
response curves are not uniform. Curves 5, 6, 7 and 9 show elements of softening 
and hardening effects. Notice that there is the same ordering of the softening 
in curves 7 and 9, with almost no softening in curve 8. The value of 
G = Cl(0)-I(Z/a)3/[II(0)-I(e^)II in eq. (6.14) for the tensile curves is shown in 
fig. 3. The effect of this term becomes important at about the yield strain. 
Thus, it is responsible for describing the majority of the macroscopic plastic 
flow. Further, shape of the curve is consistent with the step function approx- 
imation used in eq. (6.2) for determining the material parameters. The initial 
hardening rate term, in eq. (6.11) is responsible for characterizing the 
initial microstructure and influences the subsequent inelastic response. The 
transistion, from one term to the other, depends on strain rate and produces 
the various responses in calculated curves. This effect arose as a direct 
result of the potential function formulation. 

If the initial hardness, Z^, is lowered to near zero, the response becomes 
unstable between cycles of strain hardening and recovery as shown in fig. 4a. 
This effect has been observed in chemical lead at 29°C by Morrow and Halford 
(ref. 22) as shown in fig. 4b). They documented the effect as cyclic hardening 
and recovery. The calculated response in fig. 4a underscores the importance of 
the Zq and Z^ terms. 

The predicted and observed response to a hysteresis loop with a 25 second 
hold in compression is shown in fig. 5. The overall dimensions on the observed 
and predicted response are close including the compressive stress relaxation 
response. The shape of the predicted curve near yield on loading in compression 
and reloading in tension reflects the assumption of isotropic hardening. 

Obviously, there is a small element of kinematic hardening being observed in the 
response . 

Next, let us investigate the effect of initial load history on the sub- 
sequent creep response. This was done for the four initial strain rate histories 
and one creep stress as shown in fig. 6a. The predicted response curves have 

a very pronounced variation. Increasing the initial strain rate from 1.33x10 ^ 

-1 -4 

Sec to 8.33x10 in four steps produces increases in creep strain rate. 
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The observed creep response for IN100 for six levels of creep stress is 
shown in fig. 6b. It was expected that increasing stress would produce creep 
curves ordered with increasing creep rates. However, creep curve 20 is not 
ordered with respect to the other curves. If the increase creep stress of 
21 MPa between curves 19 and 20 is accompanied by a decrease in the initial 
strain rate, the observed response could be predicted. 

Further evidence of the effect of initial strain rate on the subsequent 
response is shown in fig. 7 for Rene^ 95 at 649°C. All four tests were run 
using the creep stress of 1034.3 MPa. Two tests were run on a creep frame 
with the weights placed, one at a time, on a load pan. This produced a "slow" 
initial strain rate. The other two tests were run in a hydraulic closed loop 
testing system with the initial load applied rather rapidly producing a "fast" 
initial strain rate. There is significant variation in the observed creep 
response that might be attributed to the variation in initial strain rate. 

However, the two sets of tests were run in different laboratories on different 
heats of the material. 

The effect of initial strain rate on stress relaxation was also investigated. 
Shown in Fig. 8 are four- stress relaxation curves that were calculated using 
four different initial strain rates to a strain of 0.01. The initial strain 
rate has a dramatic influence on the initial stress as might be expected from 
the tensile curves. The initial rate of stress relaxation and the final stress 
level both appear to depend on the initial strain rate. The observed stress 
relaxation from 0.01 strain is also shown. This curve is part of a creep and 
recovery history with a subsequent deformation to 0.01 strain. 

VIII SUMMARY 

The potential function formulation for the state variable evolution 
equation was based on a combined inequality that allowed for both strain hardening 
and recovery. The resulting evolution law demonstrated both of these properties. 

In particular the derived evolution equation required initial values for both the 
state variable, Z, and its rate, Z. In a given heat of material, with a fixed 
initial microstructure, the corresponding initial value for the state variable, 

Z, is assumed constant for all tests. However, the initial hardening rate, Z, 
depends on how the specimen is initially loaded. This permits’- the equation to 
account for the formation of different initial microstructures with different 
properties that can effect the subsequent macroscopic inelastic flow. 
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This effect was investigated and the results can produce response curves 
that are not similar or even ordered under various loading conditions. These 
effects were observed in IN100, Rene" 95 and chemical lead at high temperature. 
Thus, it is possible that the variations in the initial loading condition could 
account for atleast some of the variations observed high temperature tests. This 
could be particularly useful for explaining the large amount of scatter observed 
in creep tests. 
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TABLE 1. SUMMARY OF THE EXPERIMENTAL PROGRAM 


Test 

No. 

Type of 
Test 

Spec. 

No. 

Control 

Variable 

Observed 

Variable Comments 

2 

Tensile 

G1 

• -3 -1 

e «* 1.42x10 s 

o = 1116 MPa 

3 

Tensile 

T1 

• -4 -1 

X = 8.33x10 s 

Const. Hd. Rate 

4 

Tensile 

T3 

• —4 ”1 

e = 8.33x10 s 

o - 1068 MPa 

5 

Tensile 

G2 

e = 6.33xl0 _5 s _1 

a = 951 MPa 

6 

Tensile 

ENTEN 1 

e = 5.5x10 5 s 

a = 978 MPa 

7 

Tensile 

GT7 

• -5 -1 

e = 1.33x10 s 

o = 889 MPa 

8 

Tensile 

ENTEN 4 

e = 6.67x10 ^s 3 

o = 841 MPa 

9 

Tensile 

G3 

e = 1.67xl0~ 6 s -1 


10 

Creep 

Cl 

a = 496 MPa 

• -8 -1 
e = 1.8x10 s 

14 

Creep 

GT6 

o = 627 MPa 

• -8 -1 
e = 5.0x10 s 

15 

Creep 

GT5 

a = 620 MPa 

. -7 -i 

e = 1.2x10 s 

16 

Creep 

C4 

a = 696 MPa 

Relax. Obs 

18 

Creep 

GT4 

a = 827 MPa 

e = 4.5 x10~ 6 s“ 1 

19 

Creep 

C5 

a = 875 MPa 

e = 7.74xl0“ 6 s _1 

20 

Creep 

ENTEN 2 

a = 896 MPa 

e = 4.17x10 6 s 3 

21 

Combined 

C3 

e = 4.0xl0~ 3 s~ 1 

Hyst. Loop 


History 
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Figure 1 
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Figure 2. Calculated and observed tensile response curves for IN10Q at 732°C 
for several constant strain rate load histories. 
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Figure 3. Response function G = C(I(0) - I(R)3/Cl(0) - I(Eq) 3 for the tensile 
response histories shown in figure 1. 
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Figure 4. (a) Tensile responses characteristics of the model for an initial 

hardness .near zero. 

(b) Observed tensile response for chemical lead at 29°C in simple 
shear. (After Morrow and Halford L22l) . 
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Figure 5. Predicted and observed hysteresis loop response of IN100 at 732°C. 
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